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1 Introduction 



The ubiquity of linear integral equations, and in particular of Fredholm equations 
(FE) [1, 2, 3], in mathematical physics, and more broadly in analysis, cannot be 
overstated. Thus, new results in this classical field, such as the results presented in 
this paper (Eqs.(2.11) and (4.9)), should be of some interest. 

To get oriented, let us briefly recall the basic definitions and facts of Fredholm 
theory, relevant to our discussion. (We have adopted throughout this paper the 
conventions and notations of Chapter II of [3].) Thus, consider a Fredholm integral 
equation of the second type in the unknown function (j)(x), 

<P(x) = f(x) + \jN(x,y)<j ) (y)dy, (1.1) 
n 

with kernel N(x, y) and given function f(x). For simplicity, we shall take N(x, y) and 
f(x) as real functions. The generalization to the complex case is straightforward. The 
complex variable A is the spectral parameter of the equation, and Q is the domain on 
which the equation is defined. To be concrete, we shall take Q as a compact domain 
of TV-dimensional Euclidean space of volume V 

J dx = V. (1.2) 

n 

We further assume that N(x, y) is bounded on £1 

\N(x,y)\<M, (1.3) 

and that it is integrable in Q with respect to both x and y. The given function f(x) 
is assumed integrable as well. 

It is also useful to introduce the operator N and the vectors \<f>) and |/), which 
correspond to the kernel N(x,y) and functions <f>(x) and f(x). Thus, in obvious 
notations, 

N(x,y) = {x\N\y),f(x) = {x\f), and <j>(x) = {x\<j>) . (1.4) 
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In terms of (1.4), we can write Fredholm's equation (1.1) as 



(l-A#)|0) = |/). 



(1.5) 



Next, we define the n x n determinant 



( Xi x 2 ... x n \ 



N 

\Vl V2 ■■■ Vn ) 



det N(xi,yj) , 



(1.6) 



where xi, . . . y n is a set of 2n points in f2. We shall refer to the Xi as the row indices, 
and to the yj as the column indices of the symbol on the left-hand side of (1.6). In 
some of the mathematical literature (1.6) is known as the Fredholm determinant, 
however, we shall reserve this name, as is customary in vast portions of the physics 
and mathematics literature, to Fredholm's first series D(X) defined below in (1.9). 

Given (1.3), it follows from a theorem due to Hadamard that (1.6) is bounded 
according to 



By definition, D n is completely antisymmetric in the x^, and also in the y^. In view of 
(1.2) and (1.7), it is easy to see that the series (1.8) converges absolutely to an entire 
function of A. 

Fredholm's first series 




(1.7) 



Fredholm's nth minor is defined by the series 





(1.8) 




(1.9) 
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and second series 




I X S\ s 2 ... s p \ 



D(x,y;\)=N(x,y) + Y^ 



dsi . . .ds, 



V y si s 2 



s p J 



(1.10) 



correspond to setting n — and n = 1 in (1.8), respectively. 

Fredholm's first series (1.9) is, by construction, the functional determinant 



of the operator on the left-hand side of (1.5). It is usually known in the literature as 
the Fredholm determinant associated with (1.1), and we shall adhere to this convention 
here. 

From the definitions (1.8) and (1.9) we can prove the important relation 



in a straightforward manner. 

The motivation for introducing the minors stems from their important roles in 
solving Fredholm's equation (1.1) in the most general case. This is briefly reviewed 
in the appendix. In particular, for values of A such that D(\) ^ 0, the solution of 
(1.1) is determined by the resolvent kernel R(x,y;X) according to (A.l). R(x,y;X), 
in turn, is given in (A. 2) as D ^'ffl ■ For values of A such that D(\) = 0, the solution 
is given by (A. 6) and (A. 10), and involves the higher minors. 

The rest of the paper is organized as follows. In the next section we derive our new 
recursion relations (Eqs.(2.7) and (2.8)) for the minors (1.8). We then solve it and 
obtain the representation (2.11) for (1.8) as an n x n determinant over resolvents. 
In section 3 we provide an interpretation of this representation in terms of non- 
interacting fermions. Finally, in section 4 we provide an explicit formula (Eq.(4.9)) 
for the functional derivative of the nth minor with respect to the kernel. 




(1.11) 




(1.12) 
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2 The Minor D n as an n x n Determinant 



We start by deriving a couple of integral equations satisfied by D n [3]. To obtain the 
first equation, expand each of the determinants in the series (1.8) with respect to the 
row Xi, and integrate with respect to all the s- variables which occur in that term. It 
is easy to see, by elementary permutations of rows and columns, that in the pth term, 
all columns s\, . . . , s p yield the same integrated contribution. Then, after resumming 
over p, one obtains 



( %l 



V yi 



y n 



A 



£(-!)*+* jV(^,y fe )A>-i 
fc=i 



+A f N(xi,s) D n 



( X\ 



\ yi 




a/i ... 

• • • y/k 



X r 



y n 



A 



X;- 



A ds 



(2.1) 



where the symbol ^ in the upper row of -D„_i indicates that the row index Xi is to be 
omitted from the string x±, . . . , x n (and similarly for y/ k in the lower row there), and 
(Vi)s indicates that Xi in the upper row of D n under the integral should be replaced 
by the integration variable s. 

The second integral equation satisfied by D n is obtained similarly, by expanding 
each of the determinants in the series (1.8) with respect to the column yi. One obtains 

( X\ ... 




k=i 



x/ k ... 

... 4, 



X r. 



Vn 



A 



+A J N(s, yi )D n 



X r 



A ds 



(2.2) 
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We now proceed to derive our own results. We assume henceforth that D(X) ^ 0. In 
this case, according to (A.l), Fredholm's equation (1.1) has a unique solution 



<j>(x) = f(x) + \ J R{x,y-,\)f(y) dy. 
n 



(2.3) 



From this solution we construct the quantity 



/ X\ 

V yi 



X Tl 



ill ■■■ Vn 



A = A / D 



n 



( Xl 

V yi 



X r 



JJn 



A <j>{yi) dyi 



(2.4) 

Then, expand the D n under the integral on the right-hand side of (2.4) according to 

(2.2). By exploiting the fact that Xj N(x, y)(f>(y) dy = <f>(x)-f(x) = Xj R{x,y; A) f(y) dy 

a _ n 

from (1.1) and (2.3), we note that there appears a S n on the right-hand side of the 

equation, which cancels the original one on the left, leaving us with the identity 

^ X\ x n 



D r . 



Vn 



V yi ■ 

J2(-iy+ k J R{x k) y %] X) D n 



k=i 



f(yi) dyi = 
( xi ... x/ k 

\yi j/i ■■■ 



X j_ 



A f(yi) dyi 



(2.5) 



Since this identity holds for all admissible given functions f(x), we conclude that D 
must satisfy the recursion relation 

n ( X\ . . . xjk X r , 



D 



^ X ]_ ... X y t 

\ 2/1 • • • Vn 

or equivalently, 



A | ='£(-iy +k R(x k ,y i ;\)D n _ 1 
k=i 



X 

\ Vl iJi ■■■ Vn 

(2.6) 





( X\ . 


x n 


\ 


n 


( X\ 


D(X)D n 






X 


= '£(-l) i+k D(x k ,y i ;\) D n _, 






Ui • 


■ Vn 


) 


k=l 


\ yi 



Xi ... X, 



"Jk X n 



l/i ■■■ Vr, 

(2.7) 
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from (A. 2). Similarly, from (2.1), and by exploiting the associated (or transposed) 
Fredholm equation ip(x) = g(x) + A / ip(y)N(y, x) dy , we obtain the transposed iden- 
tity 



D{\)D n 



( X\ . . . x n 



V vi 



y n 



\ 



X 



= J2(-iy +k D(x z ,y k ;X) D n _, 



k=l 



( Xi ... 3?i 



A 



The two identities (2.7) and (2.8) strongly suggest that the quantity 



/ X\ 



Xr. 



V yi 



■IJn 



\ 



A 



-D r , 



J 



( X\ . 
\Vi ■ 



y n 



\ 



A 



y/k ••• y n 

(2.8) 



(2.9) 



»n-l 



A 



D(X) 

is simply the n x n determinant with entries R(xi,yj] A) and corresponding minor 

X\ . . • x/fc x n 

yi iJi ••• y n 

The proof of this proposition by induction is almost trivial: This proposition is indeed 
the content of (2.7) and (2.8) for n = 2: 

/ Xi x 2 



D(X)D 2 



A 



D(x 1 ,y 1 ;X)D(x 2 ,y 2 ;X)-D(x 1 ,y 2 ;X)D(x 2 ,y 1 ;X). (2.10) 

V yi 1/2 ' 

(For n — 1 (2.7) and (2.8) yield a trivial identity.) Then, by assuming it holds for 
A n _i, we apply it to the A n _i s which appear on the right-hand sides of (2.7) and 
(2.8), and thus observe that the latter are just the expansion of an n x n determinant 
with entries R(xi,yj] A) according to the i column and zth row, respectively. The 
proposition of the induction is thus verified for A n as well. 
Thus, we have derived our first main result: 

^ X\ . . . x n 



1 



D(X) 



D r . 



X 



det R(xi,yj; A) 



(2.11) 



\yi ■■■ y n 

Note that for n—1, (2.11) coincides with (A. 2), as it should. 
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3 Interpretation of (2.11) in Terms of Non-interacting 
Fermions 



The determinantal representation (2.11) suggests, due to Wick's theorem for non- 
interacting fermions, an interpretation of D n as the correlation function of n fermions 
and n anti-fermions. 

To this end, consider the non-interacting complex Grassmann valued field tp(x), 
living on Q, with action 

S = f ^(x) [5(x -y)- XN(x, y)] dxdy . (3.1) 



n 

Its partition function [4] is given by the path integral 

Z = J Vi\) e s = Det (l - XN) = D(X) . (3.2) 

As is well-known from the annals of quantum field theory, the non-vanishing cor- 
relation functions of (3.2) are those which contain equal numbers of ifis and ^s, 
namely, 



\ip\yi)^{xi) ■ ■ ■ {y n )^{xn) ) = / V^Vi)t s ^(y^^xx) ■ ■ ■ 4)\y n )i)(x r , 



(3.3) 



The latter are determined according to Wick's theorem as 

( V f (yi)^i) • • ■^{y n )^x n ) ) = det( ^(yiMxj) > , (3.4) 

where the two-point function is 

(^(y)rf,(x)) = (x\—^- r \y). (3.5) 
1 — AiV 

Then, observe from (3.5) and (A. 3) that 

/ N(x, z){^(y)^(z) )dz=(x |- -j- \y) = R(x, y; X) . (3.6) 

n 1 — AiV 

Thus, by linearity, we obtain from (3.4) and (3.6) that 

(^^(N^ix,) ■■ -^(y^N^M) = det R(x t ,y r , X) , (3.7) 
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where (Nift)(x) = J N(x,y)ip(y) dy. Thus, we interpret the nth minor D n , according 
to (2.11) and (3.7), as the multi-fermion correlator 



A = D(X) ( ^(y 1 )(Ni;)(x 1 ) ■ ■ ■ ^(y n ) (N^)(x n ) ) 





( X\ ■ 


x n 














\ Vi ■ 


■ Vn 





J V^V^e 8 [^(y 1 )(Ni;)(x 1 )--^^y n )(N^)(x n ) 



(3.8) 



that is, D n is the 2nth moment of the Grassmann weight e s (convoluted against n 
powers of N). As such, it might be thought of as some kind of a continuum (sup- 
plementary) compound matrix associated with 1 — AiV [5]. The latter interpretation 
might be useful in studying minors of very large order, such that the 2n points X; L and 
yi become typically dense in Q. 

As a consistency check of (3.8), let us trace it over all coordinates, and see if we 
recover (1.12). Thus, 



I 



Xj 



X dxi . . . dx n 



D(X) f {^{x 1 ){Ni,){x l )---^{y n ){N^){x n ))dx l ...dx n = 



D(X) 



i)\x)(N4))(x) dx 



(3.9) 



From (3.1), (3.2) and (3.3), we see that 



1 f 4Ym. 



D(X) \ dX j 



Thus, 





( Xi . 


x n 














V Xi . 


x n 





± 

dX, 



(3.10) 



in accordance with (1.12). 
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4 The Functional Derivative of D n with Respect to 
the Kernel 

As we have discussed above, the minors D n (1.8) determine the solution of the Fred- 
holm equation (1.1). In some applications of (1.1), the kernel N(x, y) may depend on 
a set of parameters or functions, and it may be important to determine how the solu- 
tions vary with these quantities. To this end we have first to determine the functional 
derivative of the minors D n with respect to the kernel. 

For example, in a recent paper [6], we have calculated the variation of the so- 
lution of the Gelfand-Levitan-Marchenko equation with the reflection amplitude of 
scattering theory, and deduced from it the corresponding variation of the Schrodinger 
potential and wave-function. 

The minor D n is expressed in (2.11) in terms of D(X) and R(x, y; A). It is straight- 
forward to obtain the functional derivatives of these two objects with respect to 
the kernel N(x,y) directly from (1.11) and (A. 3). Thus, consider a perturbation 
7V-> N+5N. From (1.11) we see that under this variation 8D(X) = D(X)6\ogD(X) = 
—XD(X)tr ( iJ^ ^v)? from which we infer 

= -XD(X)[5(b-a) + XR(b,a;X)} 

= -XD{X)5{b-a)-X 2 D(b,a;X). (4.1) 

Under this variation we also have R — ^ R+ ^^ ^iV ^J^ . Consequently 



SR(x, y; A) 



(x\{l + XR)\a)(b\(l + XR)\y) 



5N(a, b) 

= (8(x - a) + XR(x, a; A)) {6{b - y) + XR{b, y; A)) . (4.2) 



We can then calculate 



' D(X) det R(xi, yj; A) 



SN(a,b) V iJ 

by applying (4.1) and (4.2) as necessary. The expression we obtain in this way is 
rather cumbersome, but the plethora of terms thus obtained can be organized into a 
linear combination of the minors D n and D n ±\. 
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Instead of pursuing this line of derivation, we shall now sketch the calculation of 
the functional derivative of D n directly from (1.8), by taking the derivative of this 
series term by term. Thus, consider taking the derivative of the pth term. Let us 
split the matrix whose determinant 

( X\ X2 ••• X n S n+ i S n+ 2 • • • s n+p 



N I ds n+ i . . . ds n+p (4.3) 

& \ Hi U2 ■ ■ ■ Un s n+l s n+2 ■ ■ ■ s n+p 

is being integrated in that term into four blocks, according to NN, NI, IN and II, 
where T stands for an integrated coordinate index, and 'N' for a non-integrated one. 
Let us now scan systematically through these blocks. 

When the derivative SN ^ a ^ hits the term N(x i) yj) in the NN sector, it produces 
a factor & ^^ V b) = 5(xi — a)S(yj — b) which is multiplied by the minor of N(xi,yj), 
times a sign factor (— l) t+j . The total contribution of the NN sector to the derivative 
is the sum of all these terms: 

n 

NN = J2 {-l) i+j S{xi-a)5{ yj -b) x 



Vn Sn+1 



N | I ds n+ i . . . ds n+p 




(4.4) 



Move now to the NI block. When the derivative SN ^ a ^ hits the term N(xi, s n +j), it 
produces a factor S(xi — a)5(s n+l — b) which is multiplied by the minor of N(x i7 s n+ i), 
which contains s n+ i as a row index (but not as a column index), times a sign factor 
(— iy+ n+l . Integration over s n+ i thus replaces the row index s n+ i in that minor by 
b. Now, permute the row which used to be that of s n+ i in that minor, and move it 
in between the rows corresponding to and x i+1 . This means permuting it across 
n+l — i — 1 rows and costs a sign factor (—l) n+l ~ l ~ 1 , which combines with the previous 
sign factor to (—1), independently of I. Thus, all the p columns which intersect the 
original row Xi in the NI block make the same contribution to the functional derivative, 
after integration over the remaining p — 1 variables. Finally, summing over all the 
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Xi in the NI block we obtain the total contribution by that block to the functional 
derivative as 



NI = -p^2 d( x i - a) x 



i=l 




%n S71+I • • • S-n+p— 1 



ds n+ i . . . ds n+p -i . (4.5) 



Un $n+l ■ ■ ■ s n+p—l 

Similarly, the contribution of the entire IN block to the functional derivative is 
IN = -pJ2 5( yj -b) x 



N 



I X\ 

V Vl ■■■ (]Jj)a ... Vn Sn+l 




ri+1 • • • ds n +p—i 



(4.6) 



The last block is II. Clearly, we should discuss the diagonal terms and the non- 
diagonal terms separately. When SN ^ a ^ hits the diagonal term N(s n+ i, s n +i), it pro- 
duces a factor S(s n+ i — a)5(s n+ i — b) times the corresponding diagonal minor, which 
comes with a positive sign and does not depend on s n+ i. Thus, integration over s n+ i 
simply produces a factor 5(a — b) which multiplies the remaining integral. The latter 
is the same for all diagonal terms. Thus, the overall contribution of diagonal terms 
from the II block to the functional derivative is 

„ ( %1 • ■ ■ %n ^rt+l 

Ihiag =p5{a-b) J 



N 



n 



\ yi ... y n s n+1 




Finally, when SN s ^ a ^ hits the non-diagonal term N(s n+ i, s n+ k), I ^ k, it produces a 
factor 5(s n+ i — a)5(s n+ k — b) times the minor of N(s n+ i, s n+ k), which contains s n+ k as 
a row index (but not as a column index), and also contains s n+ i as a column index 
(but not as a row index), times a sign factor 

Integration over s n+ i thus replaces the column index s n+ i in that minor by a. 
Similarly, integration over s n+ k replaces the row index s n+ k by b. Now, move the row 
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which used to be that of s n+ k right below the row x n , and the column which used 
to be that of s n+ i immediately to the right of the column y n . These permutations 
produce a sign factor (— l) fe_1+ ' _1_1 , which combines with the previous sign simply 
to (—1). The remaining integral is independent of s& and s; and yields the same 
contribution for all the p(p — 1) non-diagonal terms. Thus, their total contribution 
to the functional derivative is 

X n b S n -|_i . . . Sn+p— 2 



ttnon-diag. = ~P(P ~ l ) j N 



( 



Xi 



ds n -\.\ . . . ds 



n+p-2 



Un a s n+l 



Sn+p-2 







Gathering all contributions (4.4)-(4.8) together, multiplying their sum by ^ and 
summing over p, we finally arrive, after some rearrangement of terms, at our second 
main result: 

r I X\ . . . X Tl 



5N(a, b) 



D r . 



>Jn 



A = 



*&1 • • • %n 

]T (-l) i+j 5(x i - a) 5( Vj -&)£> n _i ' 

ij ' =1 \ Vl ?Jj ■■■ Vn 



A 



+A^ 6(xi - a) D n 



i=i 



Xi 



2/i 



. . . (o/i)b ... x n 



y n 



A 



^ X\ ... ... ... *&Ti 



+\J2S( yj -b) D, 

3=1 



i)a ... y n 



A 



-XS(a - b) D r . 



^ X ^ ... Ob y t 



V yi 



A I - A 2 D 



I X\ ... x n b 



n+l 



V yi 



Vn a 



(4.9) 
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It is gratifying that the functional derivative of D n is expressed as a relatively simple 
linear combination of D n and D n±1 . In particular, we could think of (4.9) as recur- 
sively defining D n+1 in terms of the lower minors. This is analogous to Jacobi's recur- 
sive definition of supplementary compound matrices in the finite dimensional case [5]. 
The formuala (4.9) for the functional derivative of D n should coincide, of course, with 
the expression we would obtain by taking the derivative of D(\) dety R(xi, yf, A). 

As a simple application, let us check (4.9) for n = and n — 1. For n — 0, it 
yields 

sr>(\\ 

-X5(a-b)D(X)-X 2 D(b, a; A), 



5N(a, b) 

which coincides with (4.1). Similarly, for n — 1, we obtain 

5 ~lN(PbY = 6 ( x - a ) 6 (y- b ) D ( x ) + X6 ( x - a ) D ^y^) + X6 (y- b ) D (^ a ^) 

- X5(a-b)D(x,y;X)-X 2 D 2 ^ X y \ 

Thus, from the last two equation we can derive that 

SR(x,y;X) 6 ( D(x,y;X) \ 

SN(a,b) 5N(a,b) \ D(X) J 

5(x - a) 5(y -b) + X5(x - a) R(b, y; X) + XS(y - b) R(x, a; A) 



(4.10) 



+A 2 {R{x,y;X)R(b,a;X)-A 2 ' ' 



A (4.11) 



y a 

which coincides with (4.2) due to our determinantal representation (2.11). 

Appendix: The Solution of Fredholm's Equation 



The theory of the Fredholm's equation (1.1) and its solution is summarized in Fred- 
holm's celebrated three theorems (sometimes referred to collectively as "Fredholm's 
alternatives"). The exposition in this appendix will be very telegraphic. We will 

13 



describe the content of Fredholm's theorems in a semi-quantitative way, sufficient for 
our purposes, and refer the interested reader to the cited literature for more details. 
Broadly speaking, the solution of (1.1) depends on whether D(X) ^ or not: 
Case (1): D(X) ^ 

In this case, the solution of (1.1) involves the n = 1 minor D(x,y; A). For values of 
A such that -D(A) ^ 0, the operator 1 — XN is invertible, and (1.1) (or, equivalently 
(1.5)) has a unique solution, given by 

<t>(x) = f(x) + A f R(x, y; X)f(y) dy , (A.l) 
n 

where R(x, y; A) is the resolvent kernel of (1.1), which is given by 

R(x,y,X)^^^. (A.2) 

From (1.5) we can read off the operator R, which corresponds to R(x, y; A) = (x\R\y) 
as 

H=i(-U-lU^. (A.3) 
X\1-XN J 1-XN 

Thus, from (A.2), (1.12) (at n = 1), (1.11) and (A.3) we conclude that 

r d 
trR = / R(x, x\ A) dx = -— log D(X) , (A.4) 
J dX 
n 

which shows that the poles of R as a function of A are the zeros of -D(A) . 

Case (2): D(X) = 
If A = A such that D(X ) = 0, the homogeneous equation 

<f>{x) = X J N(x,y)<j>(y)dy (A.5) 
n 

has one or more non-trivial, linearly independent solutions 3>j(x) , % = 1, . . . is, where 
v > 1. In this case, we say that Ao is an eigenvalue of N of rank is, and refer to the 
functions <&,(#) as the characteristic functions corresponding to the eigenvalue Ao- 
(This nomenclature deviates from that of linear algebra, which would refer to ^- as the 
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eigenvalue.) Any solution of the homogeneous equation (A. 5) is a linear combination 
of the characteristic functions, i.e., the characteristic functions span Ker(l — XqN). 

As it turns out, the v characteristic functions are proportional to the z/th minor 
D u . More precisely, for a fixed set of 2v points xi, . . .y„, such that 

/ x i ■ 



D„ 



Xi. 



we have 



D„ 



\ yi 

( Xl 

V yi 



Ao \^0, 



yv 

{xii)X 

yi 



X h 



yv 



Ar 



Xl 



(A.6) 



Ao 



V yi ■■■ yv 

where the symbol (rfi)x indicates that the ith row index Xj in the upper row of the 
minor in the numerator is to be replaced by the coordinate x. The functions (A.6) 
are normalized such that 

&i(x k )=6 itk , (A.7) 

as can be seen from (1.8), since D n vanishes when any two of its row (or column) 
indices coincide. 

Since the characteristic functions are expressed in terms of D u , it does not vanish 
identically. In fact, it is the minor of lowest order which does not vanish identically as 
a function of its 2u arguments at A = Ao- In fact, (A.6) is obtained by setting n — v 
(and A = Ao) in (2.1). Furthermore, it follows from (1.12) that Ao must be a zero of 
D(X) of multiplicity which is greater or equal to u, since d ®^ might still vanish at 
A = A . 

As for solving the inhomogeneous equation (1.1) at A = A , one proceeds as fol- 
lows. First, one has to consider the transposed, or associated homogeneous Fredholm 
equation 

i){x) = A f^{y)N{y,x)dy. (A.8) 
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Since its kernel is the transpose of the kernel of (1.1), it has A as an eigenvalue of 
the same rank v. Thus, there are v independent characteristic solutions *&i(x), which 
are also expressed in terms of the i4h minor D u , similarly to (A. 6), as 

I X\ . . 



Xi 



V yi 



X 



Xj, 



An 



x v 

Vv 



(A.9) 



A 



which span Ker(l — XqN t ). 

A necessary and sufficient condition for the existence of a solution of (1.1) at 
A = A is then that the given function f(x) be orthogonal to all the characteristic 

functions ^(x), i.e., that / ^i(x)f(x) = 0, i — 1 , . . . , v. 

n 

If this condition holds, the solution (which exists) is not unique, since given a 
particular solution, one can always add to it an arbitrary solution of the homogeneous 
equation (A. 5). The piece in the general solution of (1.1) which is linear in f(x) (i.e., 
a particular solution of (1.1)) is 



+i 



M x ) = f( x ) + A o / 

n 



( X X\ 


• • • %v 


V y yi 


■ ■■ y v 


( X\ 










A 


V 2/i • • • 


y v 





Ar 



f(y) dy- 



(A.10) 
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